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Abstract. We study a class of geometric Lorenz flows, introduced independently by 
Afraimovic, Bykov & Sil'nikov and by Guckenheimer & Williams, and give a verifiable 
condition for such flows to be mixing. As a consequence, we show that the classical 
Lorenz attractor is mixing. 



1. Introduction and statement of results 

1.1. The Lorenz equations. Many systems of nonlinear differential equations that were 
first studied almost 50 years ago and which were motivated mainly by problems in geo- 
physical and astrophysical fluid dynamics and dynamical meteorology [2*llHl ll5| l2"2"]. remain 
difficult to understand rigorously to the present day. In 1963, Lorenz [T7] introduced the 
following system of differential equations: 

x = 10 (y — x) 

(1.1) y = 28x — y — xz 

z = xy — |z 

Approximate numerical studies of these equations led Lorenz to emphasise the possibility 
and importance of sensitive dependence on initial conditions even in such simplified models 
of natural phenomena. A combination of results obtained over the last 25 years PUEHHE 
|^JI2H] and culminating in the work of Tucker [2111211 gives the following statement (see 
[JHlEEj for detailed surveys): 

The Lorenz equations admit a robust attractor A which supports a "physical" 
ergodic invariant probability measure v with a positive Lyapunov exponent 

Recall that the measure v is called physical, or Sinai- Bowen-Ruelle (SRB), if for Lebesgue 
almost every solution u(t) G M 3 starting close to A and all continuous functions h : M. 3 — > M, 

1 <- T 



lim - / h(u(t))dt = / hdv. 

T->oo T J J a 

In this paper we take a further step in the understanding of the statistical properties of 
the Lorenz attractor. A measure v is mixing for a flow $ t if 

u($ t (A) n B) -> v{A)v{B) 

for all measurable sets A, B, as t — > oo. We say that the Lorenz attractor is mixing if the 
SRB measure v mentioned above is mixing. We prove the following 

Theorem 1. The Lorenz attractor is mixing. 
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In fact, our result shows that the Lorenz attractor is stably mixing: sufficiently small C 1 
perturbations of the flow are mixing. 

There are relatively few explicit examples of flows that have been proved to be mixing. 
Anosov [3] showed that geodesic flows for compact manifolds of negative curvature are 
mixing, and this was generalised (T^J to include contact flows. Moreover, mixing persists 
under C 1 perturbations. For codimension one Anosov flows [20] an d f° r Anosov flows with 
a global infranil cross-section [SJ , the set of mixing flows is C 1 open, but the corresponding 
result for general Anosov flows is not known. Recently jTT] it was shown that mixing holds 
for a C 1 -open and C r -dense set of C r Anosov flows for all r > 1. (In these references, mixing 
is proved for any equilibrium measure for a Holder potential.) However, the conditions for 
stable mixing in llUj lTT] for Anosov and uniformly hyperbolic (Axiom A flows) are not 
easily verifiable. 

The Lorenz attractor is an example of a singular hyperbolic attractor 18 j (uniformly 
hyperbolic, except for a singularity due to the attractor containing an equilibrium). Some- 
what surprisingly, we show that the singular nature of the Lorenz attractor assists in the 
search for a verifiable condition for mixing. 

1.2. Geometric Lorenz attractors. Afraimovic, Bykov and Sil'nikov pQ and Gucken- 
heimer and Williams |T51[2H| introduced a geometric model that is an abstraction of the 
numerically-observed features possessed by solutions to (jl.lj) . Tucker [2E1I21] proved that 
the geometric model is valid, so the Lorenz equations define a geometric Lorenz flow. 

Accordingly, our approach to Theorem Q is to establish mixing for geometric Lorenz 
flows satisfying certain hypotheses, and then to verify from 2fill2Tl that the hypotheses are 
satisfied for the Lorenz equations. Roughly speaking, a geometric Lorenz flow is the natural 
extension of a geometric Lorenz semiflow which is itself a suspension flow built over a certain 
type of one-dimensional expanding map /. Precise definitions are given in Section [21 In 
the literature, a standard assumption is that the map / is locally eventually onto (l.e.o.), 
see Definition El We prove that this is a sufficient condition for the corresponding flow to 
be mixing. 

Theorem 2. Let f t be a geometric Lorenz flow. Suppose that the associated one- dimensional 
map f is l.e.o. Then f t is mixing (and even Bernoulli). 

The remainder of the paper is structured as follows. In Section [21 we discuss the geo- 
metric Lorenz flow and its relations to the Lorenz attractor. In particular, we indicate how 
Theorem ^ follows from Theorem [21 In Sectional we prove Theorem [21 In Section we 
discuss extensions of our main results and some related future directions. 

2. The Lorenz attractor is a geometric Lorenz flow 

In this section we collect and organise several results from the existing literature on 
the relation between the Lorenz attractor and geometric Lorenz flows. We recall first 
of all some basic relevant facts about the Lorenz equations (jl.lj) . see [25J. The origin 
is an equilibrium of saddle type with two negative (stable) and one positive (unstable) 
eigenvalues \ ss < \ s < < X u . It is also the case that \ u > \X S \. 
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Suppose that a finite number of nonresonance conditions are satisfied so that the vector 
field is smoothly linearisable in a neighbourhood of 0. In these coordinates, the flow 
near has the form (xi,x 2 ,x 3 ) \— > (e Xut xi,e Xast x 2 ,e Xat x 3 ). By a linear rescaling, we can 
suppose that the domain of linearisation of the flow includes the cube [—1, l] 3 . Define 
E = {(x 1 ,x 2 ,x 3 ) : \x 2 \ <l,x 3 = 1} and E = {(x 1 ,x 2 ,x 3 ) : Xi = ±1, \x 2 \, \x 3 \ < 1}. 
Then we define the first hit map Pq : E — > E by 

P (xi,x 2 ,l) = (e A " r °x 1 ,e Assro x 2 ,e Asro ) = (sgnx u x 2 ,x 3 ), 

where r is the "time of flight" from E to E. Solving e XuT °Xi = sgnxi for r = —(In |xi|)/A u , 
we obtain 

P {xx, x 2 , 1) = (sgnxi, |xi| /3 x 2 , \xi\ a ), 

where a = |A S |/A U G (0, 1) and (3 = |A SS |/A U > 0. Note that P : E — > E is well-defined on 
E \ W s (0) and that r has a logarithmic singularity at xi = 0. 

Tucker j2Hll2Zj proves that there exists a compact trapping region N C E such that 
the Poincare first return map P : N \ W s (0) N is well defined. We can decompose 
P = P\oP Q where Pi is the first-hit map between E and E. Note that Pi is a diffeomorphism 
where it is defined, and that the time of flight r\ for Pi is bounded. Hence the time of 
flight r = To + T\ for the full return map P is smooth except for a logarithmic singularity 
at x\ = 0. Moreover, the following crucial hyperbolicity estimate holds. 

Lemma 2.1 ([27|). The return map P admits a forward invariant cone field. In other 
words, there exists a cone C(u) inside E at each point u G N\ iy s (0) such that (dP) u C(u) 
is strictly contained inside C(Pu). 

Moreover, there exist constants c > 0, r > 1, such that for each u G iV \ W s (0), 

\\(dP n ) u v\\>CT n \\v\\, 

for every v G C(u) and n > 1. 

A consequence of Lemma 12.11 is that the return map P has an invariant stable foliation 
with C l+£ leaves for some e > 0. Let / = [—1,1]. We obtain a singular one-dimensional 
map /:/—>/ by quotienting along stable leaves. At the same time, r reduces to a singular 
map r : / — > IR + . Let J — I \ {0}. Then it is immediate that on J: 

f(i): / is C l+£ . 

f(ii): |/( n >| > cr n for all n > 1. 
f(iii): C- 1 ^- 1 < f{x) < C\x\ a ~\ 
r(i): r(x) — >• oo as x — >• 0, 

r (ii) : \r(x) — r(y)\ < C\ In |x| — In \y\ \ for all x,y > and all x,y < 0, 

where c > 0, r > 1 are the constants in Lemma [2. 11 a = |A S |/A U G (0, 1), and C > 1 is a 
constant. Note that it follows from r (ii) that r is Lebesgue integrable and is C 1 on J. 

We have described how to pass from the original three-dimensional flow defined by the 
Lorenz equations (jl.lj) to a two-dimensional Poincare map P : E — > E and then to a 
one-dimensional expanding map /:/—>/. The process is reversible: taking the natural 
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extension of / recovers the Poincare map P and taking the suspension of the P by the roof 
function r recovers the original Lorenz flow. Alternatively, these steps can be carried out 
in the opposite order to recover the Lorenz flow as the natural extension of the suspension 
semiflow of the map / by the roof function r. The latter viewpoint is the geometric 
Lorenz flow construction [TH I2§] . For completeness, the notions of suspension and natural 
extension are recalled at the end of this section. 

Definition 1. Let /:/—>/. Assume that /(0) is undefined, with /(0 + ) = —1 and 
/(0-) = +1. It is assumed that f(l) G (0,1) and f{-\) G (-1,0). (See Figured) If 
conditions f(i) f(iii) are satisfied, then / is a Lorenz-like expanding map. 




FIGURE 1. The graph of a Lorenz-like expanding map / : [—1, 1] — > [—1, 1]. 

Definition 2. A semiflow f t is called a geometric Lorenz semiflow if it is the suspension 
over a Lorenz-like expanding map /:/—>/ by a roof function r satisfying conditions r(i) 
and r (ii) . A flow f t is called a geometric Lorenz flow if it is the natural extension of a 
geometric Lorenz semiflow. 

The results of Tucker described above imply that the Lorenz equations indeed define a 
geometric Lorenz flow. Moreover, as outlined below, Tucker showed that the associated 
expanding map / : I — ► / is transitive for the Lorenz equations, thus establishing the 
existence of the Lorenz attractor. 

In the following definition, we stress that /(0) is undefined. In particular, the endpoints 
±1 have no preimages. 

Definition 3. The map /:/—>/ is locally eventually onto (l.e.o.) if for any open set 
U C J, there exists k > such that f k U contains (0, 1). 

Remark 2.2. Since /(0, 1) contains (—1,0) and vice versa, we could equally use (—1,0) 
instead of (0, 1) in the definition of l.e.o. Clearly, the l.e.o. property implies that / is 
topologically transitive. The exact formulation of l.e.o. varies in the literature, but our 
definition agrees with the one in [T2] . 

Remark 2.3. It was anticipated in |2H| that f(ii) would hold with c = 1, r > y/2. It is 
easy to check that this is a sufficient condition for / to be l.e.o. 

Surprisingly for the actual Lorenz equations it turns out that f'(x) < 1 for certain 
values of x G I. Nevertheless, it is still the case that / is l.e.o., see Tucker \27\ . 
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It can be shown that Lorenz-like expanding maps satisfying the l.e.o. condition have 
a unique ergodic probability measure /i that is equivalent to Lebesgue (see for example 
Section EJ). The suspended measure v = \i r defines an SRB measure for the geometric 
Lorenz flow. In particular, Theorem ^ follows from Theorem 121 The latter result is proved 
in Section El 

Remark 2.4. The approach in this paper establishes weak mixing for geometric Lorenz 
semiflows. It then follows from Ratner [21] that such semiflows, and the corresponding 
flows, are Bernoulli. In particular, they are mixing. 

Suspensions and natural extensions. As promised, we end this section by recalling the no- 
tions of suspension and natural extension. Let (X, //) be a probability space and / : X — > X 
a (noninvertible) measure preserving transformation. Let r : X — > IR + be an L l roof func- 
tion. Define the suspension X r = {(x,u) : x G X, u G [0, r(x)]}/ ~ where (x, r(x)) ~ 
(fx, 0). Form the suspension semiflow f t : X r — > X r given by ft{x, u) = (x, u+t) computed 
modulo identifications. An invariant probability measure is given by fi r = fi x ij J r dm 
where I is Lebesgue measure. 

Next suppose that f t : Q — > Q is a semiflow preserving a probability measure v. Form 
the natural extension (or inverse limit) Q consisting of curves {x(s)} s >o in Q satisfying 
ft(x(s)) = x(s — t) for all s > t > 0. An invertible flow f t : Q — > Q is given by f t {x(s)} s >o = 
{x(s + t)} s >o- The projection n : Q — > Q, {x(s)} s >o i— >■ x(0) defines a semiconjugacy 
between f t and f t and there is a unique /^-invariant measure !/onfi such that tt is measure 
preserving. 

3. Mixing for geometric Lorenz flows 

In this section, we prove Theorem El Namely, we show that if / is an l.e.o. expanding 
map of Lorenz-type and r is a logarithmic roof function, then the corresponding geometric 
Lorenz flow is mixing. By Remark 12 A\ it is enough to prove that the corresponding 
geometric Lorenz semiflow is weak mixing. 

Let /:/—>/ be a Lorenz-like expanding map. As shown in Diaz-Ordaz j^j, / admits 
an induced map F that is Gibbs-Markov. More precisely there is an open interval Y G I 
containing and (modulo a set of Lebesgue measure zero) there is a partition V = {uj} of 
Y consisting of intervals, and a return function R : Y — ► N constant on partition elements 
such that the induced map F(x) = f R ^ x \x) restricts to a diffeomorphism F\ u : uj — > Y for 
each partition element uj, and such that the following conditions are satisfied: 

(a) There exists A > 1 such that {DF^l > A for each uj. 

(b) For all uj, logg^ is Holder (uniformly in uj), where g is the Jacobian of the inverse 
of F\ u : u -> Y. 

(c) R is Lebesgue integrable. 

A standard Folklore Theorem in dynamics says that conditions (a)-(c) implies the existence 
of a unique ergodic invariant measure // for / : / — > I that is absolutely continuous with 
respect to Lebesgue and whose support includes Y. In addition, the construction of 
satisfies 
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(d) G" f k u for < k < R(u) for each u. 
The following Livsic regularity result is due to Bruin et al. jH]- 

Lemma 3.1. Suppose that /:/—>/ admits an induced map F : Y —>■ Y satisfying 
conditions (a) -(d) above, with absolutely continuous ergodic measure fi. Let r : I — > R + be 
a Holder roof function with a logarithmic singularity atO. Let ip : I ^ S 1 be a fi-measurable 
function satisfying 

e ir = (^oT)f 1 , a.e. 
Then ip has a version that is Holder on Y . 

Sketch proof. We provide the details required for the reader to pass between the formulation 
in jHj and the statement of the lemma. By conditions (a)-(c), the map /:/—>/ can be 
modelled by a Young tower [31 J with base Y. The appropriate measure for the tower has a 
Holder density, so by condition (d), Bruin et al. 6, Theorem 2] applies to guarantee that %fj 
is Holder on any fixed partition element of Y and hence on the whole of Y ( (HI Remark 3] ) . 
The tower metric in [Bj coincides with the Euclidean metric when restricted to Y, so ip is 
Holder in the original metric. □ 

Remark 3.2. Related Livsic regularity results can be found in ^3]- The arguments in 
IT3] are different from each other leading to distinct results, and the approach in jB] turns 
out to more convenient for our purposes. 

Now let r : I — > IR + be a Holder roof function with a logarithmic singularity at 0, and 
let ft be the corresponding geometric Lorenz semiflow with ergodic measure /i r . 

Lemma 3.3. Suppose that the geometric Lorenz semiflow is not weak mixing. Then there 
exists a constant a > and a measurable eigenfunction ip : X — ► S 1 continuous on 

U/c>o f k Y ' > suc h th a t 

(3.1) e iar = {tfjof)^- 1 a.e. 

Proof. It can be taken as the definition of weak mixing that (j>of t — e mt <j) has no measurable 
solutions : X r —> S 1 for a > 0. Suppose that is such a measurable solution. It follows 
from Fubini that there exists e > with r > e a.e. such that o f t (x,e) = e mt (j)(x,e) for 
almost every x G X. Set t = r(x) and ip(x) = (f)(x, e). Since f r ^(x, e) = (fx, e) we obtain 
that ip o f = e iar ip. Thus ip is a measurable solution to equation (J3.1|) . By Lemma 13.11 
there is a solution ip that is Holder continuous on Y. 

It is now straightforward to show that ip is continuous on |Jfc>o f k Y- Suppose that 
z = f k y where y G Y . Since /(0) is undefined, it is certainly the case that py ^ for 
j — 0, . . . , k — 1. Hence we can choose an open set U C Y containing y such that G" pU 
for < j < k — 1. In particular, there exists 7 G (0, 1) such that ip is C 7 on K and at the 
same time e %aTk is C l on [/. 

Let Zi = f k yi where yi G U for 2 = 1,2. Applying equation (J3.1)) in the form ip o f k = 
e mrk ip, we obtain 

^(ziM^)- 1 = V(yi)^(y2) -1 e iarfc(2/l) e- iarfc(2/2) 
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so by f(ii), 

iVOOVC^r 1 ! < D\ Vl - y 2 |T < D(cA Vki - ^r. 
Hence ^ is Holder on f k U and so is certainly continuous at z as required. □ 

Proof of Theorem® Suppose that f t is not weak mixing. By Lemma l3~3l there exists a > 
and a measurable eigenfunction ip : I —>■ S 1 satisfying (J3.1j) . such that ^ is continuous on 
U/t>o f k Y- Since / is l.e.o., ip is continuous on (—1,1). In particular, ip is continuous at 

/(±1)- 

Iterate equation (j3.1|) to obtain 
(3.2) e iar e iarof = {ipo f 2 )ip-\ 

We evaluate this equation along a sequence x n > with a; n — > 0. We have the limits 
^(^n) - ► "^(0), 4>(f 2 Xn) — > 4>(f an d r(fx n ) — > r(— 1). Hence, the right-hand-side of 

equation ()3.2j) and the second factor on the left-hand-side converge as n — > oo. To obtain 
a contradiction, it suffices to choose x n — > so that e mr * Xn * ) does not converge. For n 
sufficiently large, choose b n > r(e/2 n ) such that e tabn = (— l) n . Since r is continuous on 
7 — {0} and r(x) — > oo as a; — > 0, there exists i„ > such that r(x n ) = b n . As required, 
x n — > and e tar ^ Xn ' does not converge. □ 



4. Extensions and future directions 

In this paper, we have focused on geometric Lorenz attractors satisfying a locally even- 
tually onto condition. In this section, we discuss how to relax the l.e.o. condition, and 
we describe how our results extend to larger classes of singular attractors. We end by 
mentioning some future directions. 

The l.e.o. assumption occurs in much of the literature since it is a verifiable property of 
the flow that guarantees topological transitivity of /. Moreover, it holds for the Lorenz 
equations (jl.lj) themselves P7| . However, topological transitivity for / : I — ► I is not 
crucial in this paper. We could simply work with the unique ergodic component whose 
support includes the set Y introduced in Sectional Passing to the natural extension of the 
suspension over this component yields a geometric Lorenz attractor with an SRB measure 
as before. 

The proof of mixing for this attractor also relies to some extent on the l.e.o. property, 
but it is clear that the proof goes through under a much weaker condition: namely that a 
preimage of a forward image of 1 or —1 lies in Y. (That is, there exist integers k, i > such 
that f e (l) G f k Y or f e (— 1) G f k Y.) This condition is open and dense within the class of 
Lorenz-like expanding maps, and in principle, it can be verified by a finite computation. 

Next, we note that the specific structure of Lorenz-like expanding maps is not so crucial 
for the proof of mixing. Properties f(i)-f(iii) could be relaxed, or altered completely, 
provided / admits an induced map F : Y — > Y satisfying conditions (a)-(d). The existence 
of a Gibbs-Markov induced map is a very general condition, see jSlCj- It clearly holds for 
uniformly expanding maps and, we have made use of the the recent work of Diaz-Ordaz 9J 
for Lorenz-like expanding maps. 
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Condition f(ii) in the definition of the Lorenz-like expanding map corresponds to the 
eigenvalue condition X u > \X S \ which is valid for the Lorenz attractor. The resulting class 
of geometric Lorenz attractors are sometimes called expanding. There are also contracting 
geometric Lorenz attractors |2Sj that arise when X u < |A S |. It seems likely that such 
attractors can be shown to be mixing using the ideas in this paper. 

In situations where Lemma lTTl fails. it may not be possible to reduce to a one-dimensional 
expanding map. However, it is plausible that the Poincare map P could be modelled by 
an "unquotiented" Young tower as in [3U| to which the ideas in this paper might still be 
applicable. 

In this paper, we have established mixing. An interesting open question is to prove 
results on the speed of mixing. In a different direction, it would be interesting to derive 
statistical limit laws such as the central limit theorem and invariance principles for the 
Lorenz attractor. 
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